Depth-dependent critical behavior in V2H 
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Using X-ray diffuse scattering, we investigate the critical behavior of an order-disorder phase 
transition in a defective "skin-layer" of V2H. In the skin-layer, there exist walls of dislocation lines 
oriented normal to the surface. The density of dislocation lines within a wall decreases continuously 
with depth. We find that, because of this inhomogeneous distribution of defects, the transition 
effectively occurs at a depth-dependent local critical temperature. A depth-dependent scaling law 
is proposed to describe the corresponding critical ordering behavior. 

PACS numbers: 61.05.cp, 64.60.Cn, 64.60.F-, 64.60.Kw, 61.72.Dd 



I. INTRODUCTION 

Structural defects exist in almost all real crystalline 
solids. Therefore, in order to understand structural 
phase transitions, it is crucial to understand the influ- 
ence that defects can have on ordering behavior. It has 
been shown that defects, through their accompanying 
strain fields, can change the nature of phase transitions, 
including their universal critical properties^^^. Defects 
have also been shown to be responsible for the appear- 
ance of the so-called "central peak" in diffuse X-ray or 
neutron scattering and the related "two-length scale" 
phenomena observed in an number of experimental sys- 
tems^^^^^^^^^^^. These previous studies have 
implicitly assumed that defects are homogeneously dis- 
tributed in the material, at least in the region of the 
crystal being studied, even if that region is only a "skin- 
layer", i.e., a near-surface region in which the defect den- 
sity is known to be different than in the bulk. In real 
systems, however, defects are often caused by surface 
treatments. In this case, they can be inhomogeneously 
distributed, occurring mostly in a skin-layer with a den- 
sity that continuously decays into the bulk over several 
microns ^-ii^ii^. Such an inhomogeneous distribution of 
defects complicates the ordering behavior of many real 
crystals. In this letter, using diffuse X-ray scattering in 
both reflection and transmission geometries, we analyze 
the depth-dependent critical behavior of the structural 
ordering of a crystal with this type of inhomogeneous 
defect distribution. 

Divanadium hydride (V2H) is an interstitial alloy that 
undergoes a structural phase transition between an or- 
dered monoclinic phase jSi and a disordered body cen- 
tered tetragonal phase (52 as temperature is increased (for 
the phase diagram see Ref. [HI). In the crystal we study, 
defects occur almost exclusively in a skin-layer that ex- 
tends several jiia below the surface. In the skin-layer, 
there exist walls of dislocation lines oriented normal to 
the surface^^. The density of dislocation lines in the walls 



decreases with depth. As we will see, the character of 
the structural transition in the crystal can change radi- 
cally with the depth at which it takes place. In the bulk 
material it is first order—, but it becomes continuous in 
the skin-layer and has a critical temperature and critical 
properties that depend on depth. We propose a modifi- 
cation of the scaling law for the inverse correlation length 
that accounts for its depth-dependence and allows us to 
treat the scattering measurements at different depths in 
a unified framework. 



II. EXPERIMENTAL SETUP AND RESULTS 

For our experiments we used a thin plate (0.96 mm) 
of a vanadium single crystal loaded with purified hy- 
drogen, so that it had a bulk concentration ratio of 
ch/cv = 0.525 ± 0.005. For sample preparation and 
analysis see Refs. 14 and 18. We performed X-ray exper- 
iments on the crystal in both reflection and transmission 
geometry. The reflection experiments were performed 
with MoKq;! X-rays at a rotating anode source and at 
several energies on X14A at the NSLS at Brookhaven 
National Laboratory. The high energy transmission ex- 
periments were carried out with 44.1 keV X-rays at the 
undulator beamline SRI-CAT, 1-ID, at the APS at Ar- 
gonne National Laborator}*ii. In all the experiments we 
confirmed the absence of higher harmonic contamination. 
Since we earlier observed two length scales in this crys- 
tal^^ we used the different scattering geometries and en- 
ergy ranges in order to detect separately the influences 
arising from the bulk and the skin layer. In all cases the 
sample was mounted in a strain-free manner in a vac- 
uum of ~ 10~^ torr. The temperature fluctuations of the 
entire setup were less than 0.05 K at T > 443 K. 

We have four indications of a defective near-surface 
layer on our sample: (1) A hydrogen and oxygen gradi- 
ent measured by high resolution elastic recoil detection 
analysis (HERDA) in the first 150 Ai^; (2) an oxygen gra- 
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FIG. 1: Critical temperature Tg^^ versus depth 



dient measured by secondary neutral mass spectroscopy 
(SNMS) in the first 150-200 A^^; (3) the decay of the mo- 
saic spread with depth^^^^^; (4) a larger d-spacing in the 
near-surface region than in the bulk^^. HERDA shows 
that the hydrogen content increases with depth until the 
bulk concentration is reached at a depth of about 150 A. 
In our most surface sensitive experiment we have used 
9 keV X-rays and low momentum transfers to measure 
the influence of the upper 150 A on the scattering; this 
contribution appears to be about 0.2 %, as deduced from 
the fraction R of the intensity diffracted down to a depth 
d in the sample. In a symmetric scattering geometry, R 
is given by 



R 



(1) 



where /i is the linear absorption coefficient of X-rays and 
O is the Bragg anglei^. The only defects penetrating up 
to several microns in our sample are thus the walls of 
dislocations responsible for the mosaic spread. We note 
that there is a depth-dependent stress field associated 
with this decay, caused by walls of dislocation lines^^. 

We focus here on the depth dependence of the crit- 
ical behavior in the vicinity of the /3i-/32 order-disorder 
transition in the several micron thick skin-layer. The cor- 
relation length ^ = 1/ nioT T ^ T^'^ at a constant depth 
is the inverse of the half width at half maximum of the 
critical diffuse scattering (CDS) profiles, and could be 
reliably determined only through a fit involving a convo- 
lution of the measured resolution function and the CDS 
profiles. The critical temperature T^^^ (d), whose value 
depends on the depth d (see Fig. [1]), is defined here as 
the extrapolated temperature at which the full width at 
half maximum reaches 0. Note that T^^* {d) is an "inte- 
grated" quantitity that depends on the scattering from 
the entire skin-layer down to depth d. The fact that 
Tq^^ {d) depends on depth, presumably implies that there 
is a real, local critical temperature Tc {d) that also de- 
pends on depth. Note also that the concept of a local 



FIG. 2: (Color online) Inverse correlation length k versus 
reduced temperature t = t^tI^ — 1 for different depths. The 
black circles correspond to 1.6 /xm, the red squares to 13.1 /xm, 
the blue diamonds to 18.4 /xm, the pink upward triangles to 
25 /im, the green downward triangles to 34 /xm and the orange 
stars to the bulk. 



critical temperature is well-defined only within a region 
smaller than the local correlation length. In other words, 
for any particular depth d*, when the temperature ap- 
proaches Tc (<i*) the ordering process will happen in a 
layer around that depth not thicker than the local corre- 
lation length. 

The high energy transmission experiments indicated a 
strong first order phase transition in the bulk, evidenced 
by a strong drop of the (0 5/2 5/2) superstructure in- 
tensity by a factor of more than 400 at Tg"^ (see Ref.[i3). 
This is associated with a transition width of ~ 0.3 K, to- 
gether with an abrupt broadening of the intensity profile 
at the superstructure position. 

However, in the skin layer the transition is continu- 
ous and the long-range order parameter exponent P can 
be determined from the integrated Bragg intensities of 
superstructure reflections /: 



I oc^^ = -B 



2/3 



(2) 



where ^ is the Bragg- Williams order parameter, and 
B a constant. From the corrected intensities for the 
(0 5/2 5/2) and (0 7/2 7/2) superstructure reflec- 
tions, after excluding a small two-phase region, we ob- 
tained a value of j3 = 0.18=b0.02 (see Ref. |l8|), treating 
as a fit parameter as, e.g., in Ref. HO- 
Since n scales generally as = n^t^ ^ where here 
t = — 1 is the reduced temperature^, we estimated 
the correlation length exponent u from the slope of a dou- 
ble logarithmic plot of the fitted versus t (Fig. [2]); the 
values obtained were 0.48 ± 0.05 ^ ^ 0.58 ± 0.07 when 
neglecting the crossover in u and in the susceptibility ex- 
ponent 7 at larger reduced temperatures^^. The range in 
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FIG. 3: (Color online) Proportionality factor c (eq. [3]) vs. 
depth d. The solid red line is the value of c calculated accord- 
ing to eq. [8] Inset: inverse of the mosaic spread (dimension- 
less) versus depth d. 



V is essentially attributable to experimental uncertain- 
ties, e.g., the limited resolution for smaller t. Note that 
we measure essentially the same value of v and of the 
other critical exponents regardless of depth. Thus, the 
value of the critical exponents are depth-independent. 

Our values of u and 71 (0.96 ± 0.13 as reported in 
Ref. all support tricritical behavior in the skin layer 
for small t when the tricritical point is approached along 
the T-axis. Although the correlation length decreases 
with depth, our values of u and 71 (the subscript "1" 
refers to small t) compare quite well with the theoretical 
values of u = 0.5 and 71 = 1 obtained from the analysis 
of a metamagnet which yields mean field exponents^^^^^. 
They also compare well with other systems presumed to 
be tricritical, e.g., u = 0.52 ± 0.08 and 71 = 1.05 ± 0.20 
for Nd4Cl2^. 



III. ANALYSIS AND THEORETICAL MODEL 

According to the criteria defined by Krivoglaz^^, the 
composition of the sample was close enough to the tricrit- 
ical point to observe tricritical behavior. Although the 
experimental value of /3 is smaller than the theoretically 
expected value of 0.25, it is comparable with an earlier 
measurement and with other tricritical systems if the in- 
fluence of the two-phase region is neglected^^^^^. Given 
this good agreement between the values for the critical 
exponent, we believe that tricritical behavior is dominant 
for small reduced temperatures. 

Despite the fact that the critical exponent u is depth- 
independent, there is a depth dependence to the mea- 
sured behavior of the correlation length, or, equivalently, 
to the inverse correlation length k,. To account for this 
depth-dependence, we propose that the scaling law for 



n{d) =c{d)t^ . 



(3) 



where the factor c is dependent on the depth d. To de- 
termine the function c (d) we calculated the ^/-intercepts 
of the fitted hz vs. t plots shown in Fig. [2l The resulting 
measured values for c {d) are shown in Fig. [3l 

These results indicate that, for the same change in t^, 
the inverse correlation length decreases faster for smaller 
depths. On the other hand, the critical ordering being 
measured is presumably occurring near the defect lines 
due to their strain fields^^^, which cause the appearance of 
ordered regions. We can then safely assume that the den- 
sity of ordered regions at any given depth is proportional 
to the density of dislocation lines at the same depth. The 
mosaic spread gives us a measure of this density, so we 
expect c {d) to be proportional to the inverse of the mo- 
saic spread. However, it should be noted that since hz{d) 
is calculated from the half width of the CDS profiles at 
depth it is an averaged measure. In fact, the CDS 
profiles measure ordering throughout the skin-layer to 
the depth that is probed. Thus, arguably, c{d) should 
be proportional to an integral average of the inverse of 
the mosaic spread M from the surface to depth d: 



c {d) (X 



M{z) 



dz. 



(4) 



The factor Z has to account for the absorption of the 
X-rays along the path through the material, so that it is, 
in itself, an integrated quantity: 



rKd) 

Z{d)= / e-^^dx 
Jo 



(5) 



where ja is the X-ray absorption coefficient and / is the 
effective path into the material, which, for a depth d and 
a scattering angle 1^ is given by 



l{d) 



2d 



sin 13 

We can then express Z directly as 



(d) = ^(l-e-2^^/^^"^) 



(6) 



(7) 



and, absorbing ja into the proportionality, give our final 
expression for the experimentally measured c as 



c (d) oc 



Jo M{z)^^ 
g— 2/i(i/ sin 'd 



(8) 



This quantity, calculated with average values of fi and 
sini^ and using a sigmoid function fit for M, is shown as 
a line in the main part of Fig.O The proportionality con- 
stant used, dependent on the particular sample and the 
details of the experiment, had a fitted value of 0.1 m~^. 
In the inset of the same figure we show the depth depen- 
dence of the inverse of the mosaic spread. Notice also 
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FIG. 4: (Color online) Inverse correlation scaling function 
K versus scaled reduced temperature t^, showing the data 
collapse. The black circles correspond to 1.6 /xm, the red 
squares correspond to 13.1 /xm, the blue diamonds to 18.4 /xm, 
the pink upward triangles to 25 /im and the green downward 
triangles to 34 /xm. 



that in this estimate for the experimentally measured c, 
/i and s'mi^ only play a significant role for small depths, 
due to their presence in a negative exponential. 

With this treatment for c{d)^ it is possible to collapse 
the scattering data for onto a single scaling function 



IV. CONCLUSIONS 



In summary, we have observed a change in the order of 
the phase transition in a V2H crystal, which is first-order 
in the pure bulk and continuous in a skin layer that is 
several microns thick. Defects, in the form of dislocation 
lines, exist in the skin layer and are responsible for the 
change in the ordering behavior. The density of those 
defects is inhomogeneously distributed and decays con- 
tinuously with depth from the surface. Throughout the 
skin layer, near the critical temperature, the measured 
values of the critical exponents are those of a tricritical 
point. Although the values of the critical exponents de- 
scribing the scaling behavior of thermodynamic functions 
in the skin layer do not depend on depth, the coefficients 
of the power law function of tz and, presumably, of other 
thermodynamic functions, do show a depth dependence. 
Focusing on the behavior of a^, we find that its depth- 
dependent scaling coefficient depends on an integral av- 
erage of the inverse defect density. Using this fact and 
postulating a modified scaling law for the inverse corre- 
lation length, we are able to collapse the measured data 
onto a single scaling function. Thus, we are able to treat 
scattering measurements at different depths in a single 
unified framework. 



c{d) 



(9) 



as shown in Fig. [H Note that when goes to the 
collapsed curves extrapolate to 0. This corresponds to 
a vanishing inverse correlation length, or, equivalent ly, 
to a diverging correlation length and thus suggests that 
T^^* (d) is the actual critical temperature, or at least that 
it is not too far from the real critical temperature Tc {d). 
The tails of the curves, for > 0.07 correspond to the 
critical region above the two-length scale crossover^^. 

The collapse is strong evidence that the density of de- 
fects directly drives the critical behavior in the material. 
The mechanism involved is very likely to be the inter- 
action with the stress and strain fields induced by the 
defects, which is most probably also responsible for the 
crossover to different values of the critical exponents and 
for the change in the order of the transition between bulk 
and skin-layer. 
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